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ABSTRACT 

We use galaxy-galaxy lensing to study the dark matter halos surrounding a sample of Locally 
Brightest Galaxies (LBGs) selected from the Sloan Digital Sky Survey. We measure mean 
halo mass as a function of the stellar mass and colour of the central galaxy. Mock catalogues 
constructed from semi-analytic galaxy formation simulations demonstrate that most LBGs are 
the central objects of their halos, greatly reducing interpretation uncertainties due to satellite 
contributions to the lensing signal. Over the full stellar mass range, 10.3 < log[A£,./M0] < 
11.6, we find that passive central galaxies have halos that are at least twice as massive as 
those of star-forming objects of the same stellar mass. The significance of this effect exceeds 
3(7 for log[M»/M0] > 10 . 7 . Tests using the mock catalogues and on the data themselves 
clarify the effects of LBG selection and show that it cannot artificially induce a systematic 
dependence of halo mass on LBG colour. The bimodality in halo mass at hxed stellar mass 
is reproduced by the astrophysical model underlying our mock catalogue, but the sign of the 
effect is inconsistent with recent, nearly parameter-free age-matching models. The sign and 
magnitude of the effect can, however, be reproduced by halo occupation distribution models 
with a simple (few-parameter) prescription for type-dependence. 

Key words: cosmology: observations — galaxies: haloes — gravitational lensing: weak — 
galaxies: stellar content 


1 INTRODUCTION 

The relationship between the stellar mass in galaxies and the mass 
of their underlying dark matter halos is a fundamental ingredient 
in any theory of galaxy formation and evolution. This relationship 
has been inferred using multiple observational techniques, each 
with different interpretation difficulties. These include abundance 
matching (e.g., Behroozi et al. 2010; Guo et al. 2010; Moster et al. 
2010); galaxy clustering (e.g., Zehavi et al. 2005; White et al. 2011; 
Zehavi et al. 2011; Coupon et al. 2012; Parejko et al. 2013); and 
satellite kinematics (e.g., Conroy et al. 2007 ; More et al. 2011; Wo- 
jtak & Mamon 2013). Galaxy-galaxy lensing (e.g., Hoekstra et al. 
2004; Heymans et al. 2006; Mandelbaum et al. 2006; Velander et al. 
2014; Han et al. 2015; Hudson et al. 2015; Miyatake et al. 2015; 
van Uitert et al. 2015) exploits the weak shape distortions of dis¬ 
tant galaxies due to the mass in nearby foreground galaxies, and 
therefore has the unique advantage of being directly sensitive to 
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dark matter halo masses. However, as for galaxy clustering, the in¬ 
terpretation of the results can be complicated due to the different 
lensing profiles of central and satellite galaxies (the latter of which 
include a contribution from their host dark matter halo); this is also 
the case when combining lensing with other measurements (e.g., 
Cacciato et al. 2009; Leauthaud et al. 2012a; Cacciato et al. 2013; 
Coupon et al. 2015; More et al. 2015; Zu & Mandelbaum 2015b). 

Galaxies are often classified as passive or active^, where pas¬ 
sive galaxies typically have a red rest-frame color indicative of a 
paucity of recent star formation, and active galaxies have bluer 
colors indicating recent and/or ongoing star formation. They also 
exhibit morphological differences, and variations in their distribu¬ 
tions with respect to the local density (with passive galaxies more 
prevalent in high-density environments). It is natural to ask whether 
these two classes of galaxies, which apparently have rather different 


^ Throughout this work, we use “passive” vs. ’’active” to indicate galaxy 
types, and “red” vs. “blue” to indicate rest-frame color-selected samples 
(our observational proxies for passive vs. active galaxies). 
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star formation histories, also have different relationships between 
stellar and halo mass. Unfortunately, the answer to this question 
is observationally unclear, with results pointing to differences be¬ 
tween passive and active galaxies (Hoekstra et al. 2005; Mandel¬ 
baum et al. 2006; More et al. 2011; Wang & White 2012; Tinker 
et al. 2013; Velander et al. 2014; Hudson et al. 2015; Schaye et al. 
2015), but often with quite large statistical or systematic uncertain¬ 
ties that can make it difficult to draw a definitive conclusion, espe¬ 
cially at lower mass. 

The relationship between stellar and halo mass can also be 
explored using simulation and/or galaxy modeling techniques. For 
example, semi-analytic models (SAMs) applied to subhalo catalogs 
from W-body simulations (as early as Kauffmann et al. 1997; see 
Guo et al. 2011; Benson 2012 for more recent examples) and hy¬ 
drodynamic simulations (e.g., Vogelsberger et al. 2014; Khandai 
et al. 2015; Schaye et al. 2015) predict a potentially quite compli¬ 
cated and physically-based relation between the galaxies and their 
dark matter halos, including a properly self-consistent treatment 
of the time evolution of this relationship. In both cases, there are 
parameters that must be tuned to match observations. For SAMs, 
these are explicit parameters of the semi-analytic model, while for 
hydrodynamic simulations, they model subgrid physics unresolved 
by the simulation, and again must be tuned by hand. Both methods 
are usually calibrated to match a subset of the interesting proper¬ 
ties of galaxies in real data (e.g., mass functions), and then tested 
against different properties (e.g., clustering or evolution). In con¬ 
trast, the age-matching technique (Hearin & Watson 2013; Hearin 
et al. 2014) can also be applied to catalogs from an W-body simula¬ 
tion, and assumes (rather than predicting) a particular relationship 
between the masses and ages of galaxies and those of their halos in 
order to populate the halos with galaxies. 

The various theoretical frameworks mentioned above take dif¬ 
ferent approaches to the issue of galaxy type-dependence. Obser¬ 
vational guidance on the color-dependence of the relation between 
stellar and halo mass using a simple and unambiguous measure¬ 
ment with minimal modeling assumptions may therefore be used 
to distinguish between these models. In this paper, we use a dataset 
consisting of Locally Brightest Galaxies (LBGs; Planck Collabora¬ 
tion 2013 and Anderson et al. 2015) in the Sloan Digital Sky Survey 
(SDSS) as a sample of central galaxies, classify them as passive or 
active based on rest-frame colors, and use weak gravitational tens¬ 
ing to estimate their average halo mass as a function of their stellar 
masses. We also present evidence that this LBG dataset is a rel¬ 
atively fair sample of central galaxies, and use mock catalogs to 
correct for biases resulting from the selection method and our way 
of interpreting the weak lensing measurements in terms of a halo 
mass. The use of LBGs thus reduces sensitivity to the separation 
between central and satellite contributions, which is the primary 
source of model-dependence when interpreting galaxy-galaxy lens¬ 
ing measurements, and we are left with a simple and direct mea¬ 
surement of halo masses. Wang et al. (2015) carry out a thorough 
exploration of how SAM-based catalogs capture the observed lens¬ 
ing properties of these LBGs, considering multiple cosmological 
models and SAMs; in this work, we use the mock catalog that best 
reproduces the observed properties of LBGs to interpret our mea¬ 
surements. 

The outline of this paper is as follows. In Sec. 2, we describe 
the dataset used for the measurements, and in Sec. 4 we describe the 
measurements made and how they are interpreted. Sec. 3 includes a 
description of the mock catalogs used to validate the measurement 
method, and the result of that validation. We present our results and 


compare with previous work in Sec. 5, and discuss the implications 
for galaxy formation and evolution in Sec. 6. 

All measurements made in this paper use a fiat ACDM model 
with = 0.315 in accordance with the 2013 Planck results in 
Planck Collaboration (2014), and stellar mass estimates are calcu¬ 
lated in Mq using h = 0.673 as in that work. This choice was 
made to conform with the cosmological parameters in the mock 
galaxy catalog used to validate our measurement procedure. 


2 DATA 

All data used in this paper came from the Sloan Digital Sky Survey 
(SDSS) I/II. The SDSS I (York et al. 2000) and II surveys imaged 
roughly n steradians of the sky, and followed up approximately 
one million of the detected objects spectroscopically (Eisenstein 
et al. 2001; Richards et al. 2002; Strauss et al. 2002). The imag¬ 
ing was carried out by drift-scanning the sky in photometric condi¬ 
tions (Hogg et al. 2001; Ivezic et al. 2004), in five bands (ugriz) 
(Fukugita et al. 1996; Smith et al. 2002) using a specially-designed 
wide-held camera (Gunn et al. 1998). These imaging data were 
used to create the catalogues that we use in this paper. All of the 
data were processed by completely automated pipelines that detect 
and measure photometric properties of objects, and astrometrically 
calibrate the data (Lupton et al. 2001; Pier et al. 2003; Tucker et al. 
2006). The SDSS-I/II imaging surveys were completed with a sev¬ 
enth data release (Abazajian et al. 2009), from which the LBG sam¬ 
ple described in Sec. 2.1 originates. The source catalog described in 
Sec. 2.2 uses an improved data reduction pipeline that was part of 
the eighth data release, from SDSS-III (Aihara et al. 2011); and an 
improved photometric calibration (‘ubercalibration’, Padmanabhan 
et al. 2008). 


2.1 LBG sample 

The LBG sample used in this work was initially dehned by Planck 
Collaboration (2013) and Anderson et al. (2015), so we closely 
follow their methodology with only minor modihcations explained 
below. They selected the LBGs from the flux-limited SDSS Main 
galaxy sample, using the Value Added Galaxy Catalog (VAGC; 
Blanton et al. 2005) allO sample with a constant flux limit of 
r < 17.7, including absolute magnitudes and stellar mass estimates 
using KCORRECT (Blanton & Roweis 2007), which fits stellar pop¬ 
ulation synthesis models to the five-band photometry assuming a 
Chabrier (2003) stellar initial mass function (IMF). This catalog 
was based on the seventh data release of the SDSS (Abazajian et al. 
2009). 

LBGs were selected in Planck Collaboration (2013) and An¬ 
derson et al. (2015) by defining cylinders of radius 1 Mpc, ex¬ 
tending ±1000 km/s in the redshift direction, and requiring that 
each LBG be the brightest galaxy in its cylinder according to the 
r-band absolute magnitude. In addition, a photometric redshift cat¬ 
alog (Cunha et al. 2009) was used to eliminate LBG candidates 
with a brighter companion without a spectroscopic redshift due to 
fiber collisions. The only difference between these earlier works 
and this one is that the LBG isolation criteria were previously im¬ 
posed using a WMAP7 cosmology (Komatsu et al. 2011), but here 
we use the cosmology of Planck Collaboration (2014). Imposition 
of these selection criteria left 279 343 LBGs, of which we use the 
249 818 that lie within the area covered by the SDSS source catalog 
described in Sec. 2.2. 
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We divide the galaxy sample into red and blue LBGs based on 
the g — r color K-corrected to 2 = 0.1, using 

= ( 1 ) 

for our dividing line. As discussed in Zu & Mandelbaum (2015a), 
the adoption of a stellar mass-dependent color cut to better iso¬ 
late the red sequence modifies the fractions of galaxies in the red 
and blue subsamples by at most 3 per cent for the highest and 
lowest stellar mass subsamples, so the choice of a constant divid¬ 
ing line has little impact on our results. Finally, we require red- 
shift 2 ^ 0.03, excluding the very small fraction of the sample at 
low redshift where the high galaxy flux at fixed luminosity induces 
more significant systematic uncertainty in sky subtraction and mea¬ 
surements of galaxy properties. 

Since the interpretation of stacked lensing signals can be com¬ 
plicated in cases where the mass distribution is extremely broad, 
we define relatively narrow bins in stellar mass. The halo mass dis¬ 
tribution will still be broad due to intrinsic scatter between stellar 
and halo mass, but it improves the situation compared to that with 
very broad stellar mass bins. However, relatively few high stellar 
mass LBGs are classified as blue, so we also present results for one 
very broad bin of blue galaxies above a stellar mass threshold of 
10^^ Mq. The properties of the LBG samples in each stellar mass 
and color bin are summarized in Table 1 . In addition, the number of 
galaxies and the LBG fraction as a function of redshift 2 and stellar 
mass M* are shown in Fig. 1 . The table and figure clearly indicate 
the scarcity of blue LBGs at high M*. As shown, the fraction of 
Main sample galaxies that are LBGs is highest at high stellar mass. 

Fig. 1 illustrates that the sample is quite far from volume- 
limited. We investigated the importance of the non-volume-limited 
nature of the sample by calculating LBG lensing signals for 2 < 
0.08, resulting in substantial reduction in sample size for all stellar 
mass bins except the first two (which are still not quite volume- 
limited after making this cut). However, while the errorbars in¬ 
creased significantly, there is no strong modification in the lens¬ 
ing signal, just a slight shallowing below 100 h“^kpc. For a direct 
comparison of the lensing signals between the full LBG sample and 
the 2 < 0.08 subsample, see Appendix A. As discussed in Sec. 5.3, 
these differences do not systematically change mass estimates in a 
way that suggests that a lack of volume-limited sample is driving 
our results, so for the rest of this paper we proceed with the flux- 
limited samples described above. 

Since many previous papers that have explored the relation¬ 
ship between stellar and halo mass used the MPA/JHU stellar mass 
estimates^ (an updated version of the Kauffmann et al. 2003 cata¬ 
log) instead of those from the VAGC, Table 1 includes the effec¬ 
tive mean stellar mass for each sample using the MPA/JHU masses 
(with a color- and stellar mass-dependent conversion applied based 
on comparison between the catalogs). The differences are typically 
0.1 dex for red galaxies, and slightly less for blue galaxies, with 
the correction being larger at low stellar mass. While we adopt the 
VAGC stellar masses as our canonical ones used for the analysis, 
to facilitate the comparison with previous work we will plot the 
results in Sec. 5.4 as a function of the MPA/JHU masses, denoted 


^ http : //www . mpa-garching . mpg . de/SDSS/ 


2.2 Source sample 

To measure the galaxy-galaxy lensing signal, we use a catalog 
(Reyes et al. 2012) of 1.2 background galaxies per arcmin^ with 
galaxy shape measurements (corrected for the effects of the point- 
spread function) estimated using the re-Gaussianization method 
(Hirata & Seljak 2003) and photometric redshifts from Zurich Ex- 
tragalactic Bayesian Redshift Analyzer (ZEBRA, Feldmann et al. 
2006). The catalog is characterized in detail in several papers (see 
Reyes et al. 2012; Mandelbaum et al. 2012; Nakajima et al. 2012; 
Mandelbaum et al. 2013), with well-understood systematic uncer¬ 
tainties. 


3 MOCK CATALOGS 

We use one particular Munich semi-analytic model (SAM), a 
simulation of the development of the galaxy distribution self- 
consistently including evolution and the gravitational environment, 
to create mock catalogs that are used to calibrate the mass biases 
and satellite contamination. The formation and evolution of galax¬ 
ies are modelled following the physics in the Guo et al. (2011) 
model. Central galaxies are located at the center of the friends-of- 
friends (EOF) halos, while satellite galaxies are in subhalos that 
are not at the center of FOF halos and may include orphan galax¬ 
ies for which the dark matter subhalo was disrupted before the 
present time. The stellar mass and photometric properties of model 
galaxies are determined using stellar population synthesis mod¬ 
els from Bruzual & Chariot (2003) with the Chabrier (2003) ini¬ 
tial mass function. The uncertain star formation and feedback effi¬ 
ciencies were tuned to produce close fits to the stellar mass func¬ 
tion of low-redshift galaxies from SDSS. We further used an abun¬ 
dance matching-based method to correct the stellar masses by small 
amounts in order to bring the model stellar mass function into exact 
agreement with that in the SDSS (see Wang et al. 2015 where this 
particular model is denoted Gll-P’). The model also reproduces 
the observed luminosity and auto-correlation functions of galaxies 
(Wang et al. 2015). 

The original Guo et al. (201 1) model is based on the Millen¬ 
nium Simulation (Springel et al. 2005, MS), with WMAPl ACDM 
cosmological parameters (Spergel et al. 2003). In our analysis, we 
scale the original MS to the first-year Planck cosmology (Planck 
Collaboration 2014, Ho = 67.3, = 0.315 and Ha = 0.685), 

by adopting the scaling algorithm developed by Angulo & White 
(2010) and improved by Angulo & Hilbert (2015). Thus, in the 
end, the model has the same physics and parameters as Guo et al. 
(201 1) but is based on the Planck cosmology. As shown in Angulo 
& White (2010), the rescaling algorithm can (in most respects) ac¬ 
curately reconstruct the matter field and properties of dark matter 
halos as they would have been simulated with another cosmolog¬ 
ical model; however, the concentrations of dark matter halos are 
less well-reproduced. We show the sensitivity of our results to the 
assumptions in the halo concentration vs. mass relation in Sec. 5.3. 

3.1 Mock LBG samples 

We construct mock LBG catalogs using exactly the same selec¬ 
tion criteria as for SDSS LBGs. We first project the simulation 
box along the a-axis (representing the line-of-sight direction). Each 
galaxy can be assigned a redshift based on its a coordinate and ve¬ 
locity. Selections can then be made based on the projected separa¬ 
tion and redshift difference in the same way as in the observational 
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Table 1. Summary of the properties of the galaxy samples included in each color and stellar mass bin in this analysis. The quantities that are tabulated are the 
minimum and maximum stellar mass, the number of galaxies A^gai, the effective weighted stellar mass of the galaxies taking into account their effective weight 
in the lensing measurement for our canonical (VAGC) stellar masses and the MPA/JHU stellar masses, the effective weighted redshift 2eff of the galaxies in 
the bin including the lensing weight, and the fraction of Main sample galaxies /lbg iri this bin that are selected as LBGs. 


loglo(^^) ^gal logio(^) logio(^) /lBG 


10.0 

10.4 

4244 

Red 

10.28 

10.39 

0.064 

0.13 

10.4 

10.7 

17542 


10.58 

10.70 

0.081 

0.27 

10.7 

11.0 

44724 


10.86 

10.97 

0.105 

0.47 

11.0 

11.2 

37987 


11.10 

11.20 

0.131 

0.66 

11.2 

11.4 

28008 


11.29 

11.38 

0.159 

0.79 

11.4 

11.6 

12599 


11.48 

11.56 

0.191 

0.88 

11.6 

15.0 

3195 


11.68 

11.75 

0.230 

0.91 

10.0 

10.4 

20690 

Blue 

10.24 

10.29 

0.079 

0.32 

10.4 

10.7 

30842 


10.56 

10.63 

0.100 

0.48 

10.7 

11.0 

33621 


10.85 

10.94 

0.124 

0.65 

11.0 

11.2 

11040 


11.10 

11.18 

0.155 

0.79 

11.2 

11.4 

2626 


11.28 

11.35 

0.183 

0.87 

11.4 

11.6 

325 


11.47 

11.54 

0.220 

0.90 

11.6 

15.0 

96 


11.68 

11.69 

0.246 

0.96 


data. The direct projection of the simulation box maximises the sta¬ 
tistical signal, but fails to incorporate observational effects such as 
the flux limit of the real survey, the K-corrections to obtain rest- 
frame magnitudes, and the incompleteness of close pairs caused 
by fibre collisions and the complex geometry of SDSS. Wang & 
White (2012) and Wang et al. (2014) have compared satellite prop¬ 
erties based on such direct projections of the simulation box to a 
full light-cone mock catalogue, and they found that the direct pro¬ 
jection gives unbiased results. For our later analysis we will use 
both LBGs in the z = 0 snapshot and those in earlier snapshots 
to account properly for the evolution of galaxies and their halos. 
Flowever, we have confirmed that the halo mass distributions as a 
function of stellar mass and color do not differ very strongly be¬ 
tween the z = 0 snapshot and the multi-snapshot catalogs. 

To further validate our use of these mock catalogs for this 
study, in Sec. 5.1, we will compare the lensing profiles of the simu¬ 
lated LBG samples in stellar mass and color subsamples with those 
of real LBGs. In Wang et al. (2015), the predicted dark matter pro¬ 
files from eight different semi-analytic models with varying cos¬ 
mology and model physics are compared with SDSS lensing pro¬ 
files, and this particular model is the one that gives the best match 
to SDSS. However, that comparison uses only stellar mass binning 
of the sample, so our test here using color divisions may be consid¬ 
ered a more stringent test. 


3.2 Completeness and purity of central galaxy samples 

3.2.1 Completeness 

First, we use the z = 0 snapshot of the mock catalog to assess the 
completeness of the sample of central galaxies selected using the 
LBG cut. The completeness is quantified as the fraction of central 
galaxies of a given stellar mass and color that pass the LBG selec¬ 
tion criteria. We also check whether the LBG selection is indepen¬ 
dent of the host halo mass, i.e., is the typical halo mass of central 
LBGs the same as the typical halo mass of central galaxies overall? 
Incompleteness that causes a bias in the halo mass distribution is 


more problematic for our purposes. For “typical” masses, we con¬ 
sider a straight average, a median, and since the mass 

from an NFW fit scales like on small scales (e.g., Mandel¬ 

baum et al. 2010), averaging AE corresponds to averaging 
However, since most methods measure {M 2 oom}, we will use that 
as our canonical statistic for comparison of masses throughout this 
work, and apply corrections to best-fitting masses to get an average 
halo mass. 

The completeness is shown in Fig. 2, and the ratio of typi¬ 
cal central LBG halo mass to typical central galaxy halo mass is 
shown in Fig. 3. As Fig. 2 illustrates, the completeness of the cen¬ 
tral galaxy sample selected using the LBG criterion is quite high 
(> 90 per cent) for stellar masses above 10^^ Mq, but drops be¬ 
low that. There are two reasons for this incompleteness. First, as 
discussed in Planck Collaboration (2013), many of the “missing” 
central galaxies were excluded because they are not the brightest 
object in their own halo (this explanation is also qualitatively con¬ 
sistent with the results of Skibba et al. 2011 and Hoshino et al. 
2015). Second, the typical halo masses of the central galaxies with 
stellar mass below 10^^ Mq is such that the 1 Mpc cylinder used 
in the LBG selection exceeds the virial radius. Hence low stellar 
mass centrals can be excluded because of a brighter galaxy within 
a distinct neighboring halo. This effect is more important for red 
than for blue LBGs, presumably since red galaxies tend to occupy 
more overdense regions. 

However, for this study we are more concerned with the re¬ 
sults in Fig. 3, which indicate whether the mean halo mass of cen¬ 
tral LBGs differs significantly from the mean halo mass of all cen¬ 
trals. The former is what we measure, but the latter is what we wish 
to infer. As shown, for blue galaxies the ratio of these mean halo 
masses always exceeds 0.9. For red galaxies, which suffer from 
more incompleteness, the ratio can be as low as 0.85, meaning that 
red central LBGs have masses that are biased low by about 15 per 
cent relative to those of all red centrals at the same stellar mass. 
The curves for “red” and “all” have similar trends that arise due to 
their having broader halo mass distributions than the “blue” sample, 
so even a slight mass-dependence in the LBG selection can more 
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Figure 1. Top and bottom rows show the log of the number of LBGs and the fraction of Main sample galaxies that are selected as LBGs, respectively, as 
a function of redshift and stellar mass. Results for red and blue galaxies are shown separately in the left and right columns. The color scales are defined 
self-consistently across the two columns. Horizontal lines indicate the edges of the bins in stellar mass used for our analysis. 


easily modify the mean halo mass. Moreover, the “red” and “all” 
curves are similar because the contribution of red LBGs to the “all” 
curve is boosted by the overall higher halo masses for red LBGs. 
If the ratios in Fig. 3 had been very far from 1 (e.g., 0.5) then this 
would call our entire method into question; but we will correct for 
this bias using the curves in Fig. 3, relying on the realism of the 
mock catalogs to provide valid corrections. 

We emphasize that Fig. 3 rules out the idea that we might ob¬ 
serve different halo masses for the red vs. blue LBG samples solely 
due to selection effects (e.g., that the red LBG sample might be 
heavily biased towards high-mass halos due to the selection pro¬ 
cedure preferentially excluding low-mass red centrals that are in 
overdense regions). While the corrections derived from Fig. 3 are 
in principle dependent on the details of the mock catalogs, we con¬ 
firm in Sec. 5.3 that even significant modifications in the mock cat¬ 
alogs do not lead to major changes in these corrections; alternative 
results for two different mock catalogs discussed in that section are 
shown on Fig. 3. 

Our finding that there is a small bias in the mean halo mass 
due to LBG selection is consistent with the findings of More 
et al. (2009). They compare the results of an iterative, luminosity- 


dependent isolation criterion with a fixed aperture criterion (which 
is stricter than the one used here, as it requires a factor of two dif¬ 
ference in luminosity between the LBG and anything else in its 
cylinder, instead of that it be strictly the brightest). As shown in 
their figure Al, the fixed isolation criterion can give a stellar mass- 
dependent suppression of the mean halo mass of selected galaxies 
compared to that of all central galaxies at that stellar mass. While 
a quantitative comparison is not possible due to the different cylin¬ 
der parameters and brightness requirements, our results are quali¬ 
tatively in agreement with theirs. 

3.2.2 Purity 

In addition, it is important to understand the purity of the LBG 
sample as a set of central galaxies. Contamination from satellites 
can modify the shape and increase the amplitude of the weak tens¬ 
ing signal on scales above ~ 300/i“^kpc (e.g., Mandelbaum et al. 
2005b), affecting mass estimates. This effect is particularly impor¬ 
tant if the contaminating satellites reside in very high-mass halos. 
Thus, Fig. 4 shows the purity of the LBG sample in the z = 0 
snapshot of the simulation, defined as the ratio of number of central 
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Figure 2. The completeness of the central galaxy selection as carried out 
using the LBG criterion in stellar mass bins, for red, blue, and all LBGs 
shown as red, blue, and black points, respectively. 


LBGs to number of LBGs as a function of stellar mass and color. 
As shown, at stellar masses below 10^^ Mq, the purity of the cen¬ 
tral sample determined using the LBG selection exceeds ~ 90 per 
cent. For higher stellar mass, it can go as low as 82 per cent. This is 
most likely related to the fact that (as mentioned in Sec. 3.2.1) the 
brightest galaxy in a halo is not always a central galaxy. 

However, the magnitude of the impurity alone cannot he used 
to quantify how the contaminating satellites bias the mass estimates 
from fits to AE. The bias in the masses depends on the details of 
the subhalos and larger host halos in which the contaminating satel¬ 
lites reside, the distribution of their separations from the center of 
the host halo, and other effects. For example, in the case of con¬ 
tamination due to the brightest galaxy not being the central, it may 
often be the case that the mass of the host halo is not very differ¬ 
ent from the mass one would expect if it had been a central, so the 
contamination in the lensing signal and mass estimates will not be 
very large. In the next section, we will quantify the implications 
of satellite contamination on mass estimates directly, by fitting the 
lensing signal in the mock LBG samples (including the contami¬ 
nating satellites) to estimate masses, and directly determining the 
bias. 


4 METHOD 
4.1 Measurements 

Galaxy-galaxy weak lensing, the coherent tangential shape distor¬ 
tion (or shear) of background galaxies due to the matter in fore¬ 
ground lens galaxies, probes the connection between the lens galax¬ 
ies and matter via their cross-correlation function £,gm- This cross¬ 
correlation can be related to the projected surface density^ via 

=pj [l + (y^rl q-TT^^] drr, (2) 


^ Here we neglect the very broad radial window function, given that it 
makes sub-0.2% differences in the signals on the scales used in this work, 
and neglecting it enables faster computations of theoretical signals. 


Figure 3. The ratio of the mean halo mass of central LBGs to the mean 
halo mass of all centrals in stellar mass bins, for red, blue, and all LBGs 
shown as red, blue, and black points, respectively. The additional magenta 
and cyan points labeled as “alt” in the legend were calculated using alternate 
mock catalogs that will be discussed in more detail in Sec. 5.3, as a way of 
quantifying systematic uncertainty in these results. 



Figure 4. The purity of the central galaxy sample obtained using the LBG 
criterion in stellar mass bins, for red, blue, and all LBGs shown as red, blue, 
and black points, respectively. 


where p is the mean matter density, is the physical projected 
separation and tt is the line-of-sight distance from the lens. The 
surface density is then related to the observable quantity for lensing, 
the surface density contrast: 

AE(rp) = 7t(rp)E„,t = S(< r^) - E(rp). (3) 

Here E(< rp) is the average surface density within rp. The observ¬ 
able AE is the product of two factors, a tangential shear yt and a 
geometric factor 

^ 4^ DlDls 

where Dl and Ds are angular diameter distances to the lens and 
source, Dls is the angular diameter distance between the lens and 
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source, and physical (rather than comoving) coordinates are used 
for all quantities throughout this paper. 

The lensing measurement begins with identification of lens- 
source pairs, with source photometric redshift larger than the lens 
spectroscopic redshift. Inverse variance weights are assigned to 
each lens-source pair, including both shape noise and measurement 
error terms in the variance: 


Wla = 




(5) 


where (j\ is the shape measurement error due to pixel noise, and 
(j|n is the root-mean-square intrinsic ellipticity (both quantities are 
per component, rather than total; the latter is fixed to 0.365 follow¬ 
ing Reyes et al. 2012). Use of photometric redshifts which have 
nonzero bias and significant scatter results in a bias in AE that can 
be easily corrected using the method from Nakajima et al. (2012). 
This bias is a function of lens redshift, and is calculated includ¬ 
ing all weight factors for each lens sample taking into account its 
redshift distribution. 

AE in bins in can be computed via a summation over lens- 
source pairs “Zs” and random lens-source pairs “rs”: 


AE(rp) 


2 ^ 


( 6 ) 


where et is the tangential ellipticity component of source galaxy 
with respect to the lens position, and TZ is the shear responsivity 
(Bernstein & Jarvis 2002) that converts from the ensemble average 
distortion to shear. The division by Wrs accounts for the fact 
that some of our “sources” are physically associated with the lens, 
and therefore not lensed by it (see, e.g., Mandelbaum et al. 2005a). 
Finally, we subtract off a similar signal measured around random 
points with the same area coverage and redshift distribution as the 
lenses, to subtract off any coherent systematic shear contributions 
(Mandelbaum et al. 2005a); this signal is statistically consistent 
with zero for all scales used in this work. 

To calculate the error bars on AE, we use the bootstrap resam¬ 
pling method, dividing the lens sample into 100 patches on the sky. 
Previous work (e.g., Mandelbaum et al. 2005a, 2013) has shown 
that for the scales used here, the errorbars are dominated by un¬ 
correlated shape noise (whereas on larger scales, correlated shape 
noise, systematics, and cosmic variance result in non-negligible 
correlations between bins in Vp). For this reason, we use only di¬ 
agonal errors throughout this work. The maximum scale used for 
the fits in the lensing analysis is 1 Zi“'^Mpc, which for a typical 
lens redshift is far below the typical size of each jackknife resam¬ 
pled region. Thus, the bootstrap method is a reasonable approach 
to getting the covariance matrix for the projected mass profile, and 
has been shown to agree with several other approaches in the shape 
noise-limited regime (Mandelbaum et al. 2005a). 


4.2 Mass estimation 

When estimating masses, we take advantage of the simplicity of 
the LBG sample as a nearly pure sample of central galaxies. All 
of the projected mass around these lens galaxies contributes to the 
galaxy-galaxy lensing signal. This includes contributions from the 
host dark matter halo in which the lens galaxy resides (“1-halo 
term”), and from other dark matter halos (“two-halo term”) that 
are part of large-scale structure associated with the lens. For central 
galaxies, the 1-halo term simply corresponds to the AE for the host 
dark matter halo. For satellite galaxies, there are two contributions 
to the 1-halo term: on small scales, a contribution from the satellite 


subhalo, and on larger scales (0.3-2Zi“'^Mpc) a contribution from 
the host halo itself. For a nearly pure sample of central galaxies, 
we can restrict ourselves to the 1-halo regime and take advantage 
of the fact that the halo profiles of centrals are approximately NFW 
(Navarro et al. 1997) profiles. 

For all but the highest mass bin, we use a minimum radius 
of rp = 50Zi“^kpc (physical), whereas for the highest mass bin 
where there is evidence of large-scale light profiles interfering with 
the measurements of the nearby fainter galaxies used for lensing 
measurements, our minimum scale is rp = 70Zi~^kpc. The contri¬ 
bution of the galaxy stellar mass to the lensing signal is negligible 
compared to the dark matter halo contribution for Vp above our min¬ 
imum value. In all cases, the maximum scale is 1 Zi“^Mpc, ranging 
from several times the virial radius to only a small amount outside 
the virial radius. We rely on the calibration factors derived from us¬ 
ing the same fitting procedure on mock catalogs to correct for any 
small contamination of the signal on these scales by the two-halo 
term, which may be more important at lower mass (where we go 
farther beyond the virial radius). 

Our fits have a fixed concentration-mass relation, correspond¬ 
ing to that from Diemer & Kravtsov (2015), and evaluated at the 
lensing-weighted mean redshift of each sample. However, when 
presenting our results we will estimate the sensitivity of our results 
to assumptions about the concentration-mass relation. The masses 
quoted in this paper correspond to a spherical overdensity of 200 
with respect to the mean density (M 2 oom), using Planck Collabo¬ 
ration (2014) cosmological parameters (Qm = 0.315). 

As shown in Mandelbaum et al. (2005b), using fits to NFW 
profiles when the signal is actually the average of many NFW pro¬ 
files spanning a range of masses and concentrations can give a best¬ 
fitting mass that is biased with respect to the true mean mass by a 
few to tens of per cent (depending on the width of the mass distri¬ 
bution). This point is quite relevant given that Planck Collaboration 
(2013) showed that the typical FWHM of the halo mass distribution 
at fixed stellar mass for LBGs is ~ 1 dex. In addition, NFW fitting 
can result in mass biases if the assumed concentration is incorrect, 
or if the fit includes scales outside the virial radius where the simple 
extrapolation of the NFW profile may be incorrect. Our approach in 
this work is to use A-body and SAM-based mock catalogs to cali¬ 
brate out such biases as well as those due to satellite contamination 
in our LBG sample (see Sec. 3). 

Errorbars on the mass estimates are calculated using the boot¬ 
strap. Each bootstrap-resampled AE is fit using the procedure de¬ 
scribed in this section, and the quoted results are the median of the 
mass distribution from all bootstrap-resampled datasets, with error- 
bars indicating the 16 and 84 percentiles. For all but the lowest mass 
bin (for red and blue) and the two highest-mass blue samples, where 
the signal-to-noise is lower, using these percentiles gives nearly 
identical errors to using the standard deviation. Also, the median 
of the bootstrap masses is consistent with the mass obtained by di¬ 
rectly fitting the average signal to within < O.lcr. It should be kept 
in mind that we are quoting percentiles of the distribution of mean 
halo mass values from the bootstrap; these are errors on the mean 
mass, not a measure of scatter in the halo mass distribution (which 
is much broader than our quoted errors on the mean). 

To summarize, we are inferring {M 2 oom)(AT*), the average 
halo mass in bins of stellar mass. This is the most observationally 
accessible quantity in a lensing measurement. Other works describe 
the relationship between stellar and halo mass in other ways. For 
example, halo models commonly infer a stellar-to-halo mass rela¬ 
tion (SHMR) parametrized as Mt{M 2 oom) which is a continuous 
function representing the median stellar mass of central galaxies 
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occupying halos of fixed halo mass. In principle, the quantity that 
we measure in this work for red and blue central galaxies depends 
on all of the following: the continuous SHMR, the size of the in¬ 
trinsic and observational scatter about the median SHMR, the halo 
mass function, and any halo mass dependence in the fraction of 
red vs. blue galaxies. It is important to bear in mind the distinction 
between the observationally-accessible {M 2 oom}{Mt) and the un¬ 
derlying SHMR when comparing our results with others that infer 
a parametrized M, {M 2 oom)- 


4.3 Fidelity of mass estimation procedure 

Next, we consider the fidelity of the mass estimates from our NFW 
fitting procedure. The average of a collection of NFW profiles does 
not look precisely like an NFW profile with the average mass, and 
simplifications in the assumed mass-concentration relation along 
with deviations from an NFW profile extrapolated beyond the virial 
radius can cause biases. To quantify these effects, we applied our 
mass fitting procedure to stacked AE profiles for central LBGs 
selected across multiple redshift snapshots in the mock catalogs. 

A second complication is that satellites that are accidentally 
included in the LBG sample, which are not completely negligible 
(Sec. 3.2.2), modify the AE profiles on intermediate scales and 
bias the mass estimates. To address this point, we applied our mass 
fitting procedure to stacked AE profiles for all (not just central) 
LBGs selected across multiple redshift snapshots while matching 
the joint redshift-stellar mass distribution of the real LBG sample, 
and including redshift-dependent weights consistent with the tens¬ 
ing analysis in the real data as in Wang et al. (2015). This test also 
includes the issues covered by the first test, but carrying out both 
tests separately allows us to assess the relative importance of satel¬ 
lite contamination vs. other issues. 

Results are shown in Fig. 5, where we plot the bias of the best¬ 
fitting masses with respect to the true average halo mass for central 
galaxies in this stellar mass bin, defined as the ratio of the best¬ 
fitting mass to the average central galaxy halo mass, minus 1. For 
the overall LBG sample, this is denoted Mfa /— 1, while 
for central LBGs it is denoted — 1. As shown, 

the bias tends to be negative, i.e., an underestimation of the mean 
mass. This is also true if we use the lensing signal-weighted average 
mass ((ATjoq^)^^^), which makes the bias slightly less severe. Our 
finding that the best-fitting mass is below the mean mass is consis¬ 
tent with the results of Mandelbaum et al. (2005b), who found that 
when assuming a single value of mass, the best-fitting mass tended 
to be below the mean mass and above the median mass. In that 
work, the bias was worse when the mass distribution was broader. 
Fig. 6 shows that this is the case here too, though the trend with 
breadth of the halo mass distribution is weak; the bins in Fig. 5 that 
have a larger bias for blue LBGs at high stellar mass have a broader 
halo mass distribution, also resulting in a bigger difference between 
(Maoo^) and 

The difference between the results for the two tests shown 
in Fig. 5 is very small, indicating that satellite contamination is 
not significantly biasing the mass fits, perhaps for the reasons sug¬ 
gested in Sec. 3.2.2. Moreover, the biases in the masses for red and 
blue galaxies are fairly similar, so they should not compromise our 
ability to distinguish between the halos of red and blue galaxies. 



Figure 5. The bia.s in the best-fitting NFW mass compared to the mean halo 
mass for central LBGs, for red (top), and blue (bottom) LBGs respectively. 
The two lines that are shown indicate the bias in the estimated masses of 
central LBGs, which result from the fitting procedure, and in the estimated 
masses of all LBGs (including those that are not central, thus including the 
effects of satellite contamination). The ideal unbiased case is shown as a 
horizontal dashed line. 

5 RESULTS 
5.1 Lensing signals 

Figure 7 shows the lensing signals for the LBGs in stellar mass 
bins, for red and blue LBGs'^. We first consider the basic question 
of consistency between the signals for red and blue galaxies for 
50 ^ R ^ 1000 /i“^kpc. We carry out a non-parametric test 
for consistency, defining 

2 (AErcd — AEbluc)^ 

X = Z. ^2-T^2- (7) 

i ^^red ^^blue 

where the summation is carried out over bins in Vp. This test is 
interesting when the stellar mass distribution for red and blue LBGs 
within the bin is sufficiently similar that differences in AE cannot 
come just from the different stellar mass distribution (this condition 
is violated by the broad M* > IO^^Mq bin, which is excluded 
from this test). 

Table 2 shows the degrees of freedom, and p-values for 
each stellar mass bin and combined across all stellar mass bins. 
The number of degrees of freedom appears larger than expected 
from the number of data points in Fig. 7 because the data were re¬ 
binned in the figure to make it easier to read. As shown, adopting 
the convention that p < 0.05 indicates a lack of consistency be¬ 
tween the red and blue samples, there is a clear inconsistency for 
three of the bins, and a marginal inconsistency (p = 0.06) for a 
fourth bin. In the other bins, the errorbars are sufficiently large for 
the blue LBG sample that it is difficult to draw any conclusion. The 
over all seven bins is p = 1.2 x 10““^, ruling out the null hy¬ 
pothesis that red and blue galaxies at all stellar masses in this work 
have consistent distributions of halo masses. 

As another diagnostic of the lensing profiles, we also consider 
the lensing profiles of all SDSS Main sample galaxies using the 

The data shown in Figs. 7-9 can be downloaded in machine-readable for¬ 
mat fromhttps ; //github . com/rmandelb/mandelbaum-data. 
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Figure 6. Top: The ratio of the lensing signal-weighted average mass 
to the mean mass {M 200 m) for central LBGs in the mock 
catalog. Bottom: The width of the halo mass distribution, quantified as the 
ratio of the halo masses corresponding to the 84th and 16th percentile of the 
central LBG halo mass distribution. 


same stellar mass and color divisions, but without the LBG se¬ 
lection imposed. This provides an alternative means (besides the 
mock catalogs used in Sec. 3) of checking whether the LBG se¬ 
lection chooses a biased subset of the central galaxies. Ideally, the 
characteristic “bump” that satellites can cause on ~ 300 fe“^kpc to 
1 /i“^Mpc scales should disappear when comparing the signal for 
all Main sample galaxies against that for LBGs, but the signal on 
smaller scales should be essentially preserved. This comparison is 
shown in Figs. 8 and 9 for blue and red LBGs, respectively. Fig. 8 
uses a single wide bin at high stellar mass due to the small number 
of blue galaxies above M* > 10^^ Mq. As shown, the signals for 
blue Main sample galaxies and LBGs in a given stellar mass bin 
are nearly identical, with the LBGs showing some signal deficits 
above ~ 300 /i“^kpc as expected. The differences are small pre¬ 
sumably because only a small fraction of blue galaxies tend to be 
satellites in higher mass halos (e.g., Mandelbaum et al. 2006). For 
the red galaxies, as shown in Fig. 9, the signals for Main sample 
red galaxies and LBGs at fixed stellar mass are quite consistent be¬ 
low ~ 300 /i~^kpc, and show a characteristic difference above that 
indicating that many Main sample red galaxies are satellites (espe¬ 
cially at lower stellar mass), and that the LBG selection is effec- 


Stellar mass range 

dof 


p-value 

[10,10.4] 

15 

19.7 

0.18 

[10.4,10.7] 

15 

18.8 

0.22 

[10.7,11] 

15 

26.8 

0.03 

[11,11.2] 

15 

24.2 

0.06 

[11.2,11.4] 

15 

29.8 

0.01 

[11.4,11.6] 

15 

28.7 

0.02 

> 11.6 

15 

18.8 

0.22 

All 

105 

166.9 

1.2 X 10“"^ 


Table 2. Results of a test for consistency between the lensing signal for 
red and blue LBGs in each stellar mass bin, and across all stellar masses. 
This test uses the signals from 50-1000 /r“^kpc. 

lively removing most or all of the satellites. This result provides a 
separate validation of the LBG criterion as a way of selecting a fair 
sample of central galaxies, in addition to the mock catalog-based 
validation from Sec. 3. 

Finally, as a validation of our use of mock catalogs to correct 
for biases in the fitting procedure, Fig. 10 shows, for all stellar mass 
bins and both colors, a comparison between the real LBG AE and 
the LBG AE profiles in the mock catalogs. These were computed 
via direct cross-correlation of the simulated LBGs and the (sub¬ 
sampled) dark matter density distribution in the simulations. It is 
important to keep in mind that these are not fits, and that the semi- 
analytic model in the mock catalogs was not specifically tuned to fit 
lensing profiles at all. With that caveat in mind, the agreement be¬ 
tween the data and mock catalogs is quite remarkable especially for 
the bins with M* < Mq. While the signal for blue galax¬ 

ies in the mocks is slightly high compared to that in the data, the 
match is in general close enough that the mocks cannot have too 
different a halo mass distribution (as a function of stellar mass and 
color) compared to the real data; for a more detailed exploration of 
this point, see Wang et al. (2015). The possible exception to this 
statement is in the bins with M, > Mq, in which there are 

almost no lenses in the real data, and our constraints are therefore 
extremely weak. The red LBG data in the highest-mass bins show 
more tension, at roughly the 20 per cent level (30 per cent in mass). 


5.2 Mass estimates 

In this section, we present mass estimates for the red and blue LBG 
samples. Results are shown in Fig. 11. For the two lowest mass 
samples, the mass difference (red — blue) is consistent with 0 at 
approximately the 1.5a level, while for the higher stellar mass sam¬ 
ples, the difference is positive at > 2a. 

The results in Fig. 11 include the corrections from mock cat¬ 
alogs necessary to ensure unbiased recovery of the average central 
halo masses, described in Sec. 3. In particular, we include two sep¬ 
arate corrections: a color- and stellar mass-dependent correction 
for biases in the fitting procedure in recovering (M 2 oom) (due to 
basic inadequacy of the NFW fitting, plus satellite contamination; 
Sec. 4.3), and a correction for the fact that the central LBG sample 
has a slightly different average mass compared to the full central 
galaxy sample (Sec. 3.2). As shown previously, these corrections 
go in the same direction for red and blue galaxies, with slightly 
different magnitude. The top panel of Fig. 11 also shows the re¬ 
sults without the corrections; it is clear that these corrections cause 
> la changes in the mass estimates, but are not the reason for our 
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Tp [/i^^Mpc] Tp [/i^^Mpc] Tp [/i^^Mpc] Tp [/i^^Mpc] 

Figure 7. Lensing signals for LBGs in stellar mass bins, for red and blue LBGs shown as red and blue points with en'ors, respectively. The thick lines show 
the best-fitting NFW profiles, while the vertical dashed lines show the lower and upper limits in Vp used for the fits. For viewing purposes, the data for blue 
LBGs in the four narrow bins with M* > 10^^ Mq have been rebinned into bins that are twice as broad as those shown for red LBGs, due to the lower 
signal-to-noise ratio. 



Figure 8. Lensing signals for blue Main sample galaxies vs. blue LBGs in stellar mass bins shown as black and light blue points with errors, respectively. The 
vertical dashed lines show the lower and upper limits in rp used for the fits of the LBG signals to NFW profiles. 
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Figure 9. Lensing signals for red Main sample galaxies vs. red LBGs in stellar mass bins shown as black and red points with errors, respectively. The vertical 
dashed lines show the lower and upper limits in Vp used for the fits of the LBG signals to NFW profiles. 


conclusion that red and blue central galaxies have different average 
halo masses in fixed stellar mass bins. 

The corrected values of mass and their uncertainties are tabu¬ 
lated in Appendix B, Table B1 . 

5.3 Systematic error budget 

Here we discuss the primary sources of systematic uncertainty in 
our results. As discussed in Mandelbaum et al. (2013), a reason¬ 
able estimate of the multiplicative calibration uncertainty in the 
lensing signal using this catalog is at most® 5 per cent, when con¬ 
sidering shear systematics and photometric redshift systematics to¬ 
gether. Since NFW mass fits scale like on small scales, this 

implies an 8 per cent overall mass calibration uncertainty at the 2a 
level, correlated between all the samples (i.e., it would cause co¬ 
herent shifts of all results in the same way). 

However, Fig. 11 clearly reveals a more important source of 
systematic uncertainty. As shown, the correction factor that we ap¬ 
ply to remove the effects of incompleteness, impurity, scatter in 
the mass distribution and other uncertainties in our modeling pro¬ 
cedure is significantly larger than the statistical error for nearly all 

® This should be considered more like a 2a range than a Itr range. 


the red samples, and corresponds to roughly la statistical errors for 
most of the blue samples. Thus, the validity of this correction fac¬ 
tor, which encodes our understanding of multiple different issues as 
outlined in Sec. 3, is our primary source of systematic uncertainty. 
As already noted in Sec. 5.1, there are several empirical tests of 
the lensing signal that validate our conclusions from the mock cat¬ 
alogs, so we have reason to believe the uncertainty is significantly 
smaller than the size of the correction itself. 

The corrections that were determined in Sec. 3 are a combi¬ 
nation of the mass incompleteness correction, and a correction to 
our best-fit masses to estimate the mean mass (given the halo mass 
distribution, which if broad, can result in significant differences be¬ 
tween best-fit and mean masses). The results in Wang et al. (2015) 
provide some insight into the uncertainty in the second part of the 
correction factors. In that work, multiple different mock catalogs 
were used to model the LBG lensing and clustering signals (using 
the same stellar mass bins, but without color division). The halo 
mass distribution were shifted by an overall multiplicative factor 
to match the lensing signals for LBGs in stellar mass bins, and 
the mean masses for all models (after this shift) were compared to 
get a scatter on the mean. This scatter in the mean mass reflects 
the different shapes of the halo mass distributions in the differ¬ 
ent models. For the seven stellar mass bins, excluding the broad 
one with M* > 10^^ AT©, the scatter in the base-10 logarithm of 
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rp [A^'Mpc] r, [ft^'Mpc] Tp [h^'Mpc] i-p [/i^'Mpc] 


Figure 10. For all stellar mass bins as labeled on the plot, we compare the leasing AS profiles for red and blue LBGs in the data (the points with errors, same 
as in Fig. 7) with those in the mock catalogs (the solid lines). The vertical dashed lines show the lower and upper limits in Tp used for the fits of the LBG 
signals to NFW profiles. 


the mean halo masses across all the mock catalogs is 0.035, 0.028, 
0.069, 0.055, 0.060, 0.081, and 0.11 after correcting to our adopted 
halo mass definition. Adopting these as systematic uncertainties is 
somewhat conservative, given that the extreme cases that drive the 
scatter to large values had halo mass distributions that may have 
resulted in an initial failure to reproduce LBG lensing signals. 

Table B1 shows that if we assume the same uncertainties for 
both colors, the ratio of (conservative) systematic to statistical un¬ 
certainty for red central LBG {M 2 oom) in our stellar mass bins is 
0.2, 0.2, 1.8, 1.4, 2, 2.8, and 2. For blue central LBGs, these ra¬ 
tios (using the smaller upper errorbar) are 0.2, 0.2, 0.8, 0.6, 0.4, 
0.2, 0.2. Thus we conclude that our conservatively-adopted sys¬ 
tematic uncertainty is subdominant to statistical errors in all cases 
for the blue LBGs, while for the red ones, they are comparable to 
the statistical errors above M* > 10'^°'^ Mq. However, even in the 
highly unlikely case that the results shift by the entirety of that sys¬ 
tematic uncertainty in opposite directions for red and blue galaxies, 
this would not remove the majority of the observed mass difference 
between red and blue central galaxies. 

To assess the systematic uncertainty in the incompleteness 
corrections, which account for the fact that the mean halo masses 
of central LBGs differ from those of all central galaxies, our pro¬ 
cedure is as follows. We adopted two other semi-analytic models 
from Guo et al. (201 1) and Henriques et al. (2015), which are based 
on the WMAPl (Spergel et al. 2003) and Planck (Planck Collab¬ 
oration 2014) cosmologies, respectively. The former has a similar 
semi-analytic prescription as the mock catalog used throughout this 


work, whereas the latter follows a different prescription (for more 
detailed comparison between the LBG samples in these mock cata¬ 
logs, see Wang et al. 2015). These models were specifically chosen 
due to the fact that their halo mass distributions at fixed stellar mass 
are notably different from those in our fiducial mock catalog, result¬ 
ing in a substantial mismatch when compared with the observed 
lensing signals (unlike for our fiducial mock catalog), and provid¬ 
ing a quite conservative upper bound on systematic uncertainty in 
the incompleteness correction. For these models, the completeness 
correction for blue (red) LBGs ranged from 1-4 per cent lower (5 
per cent lower to 8 per cent higher) than our fiducial one, with some 
stellar mass dependence in both cases. Fig. 3 includes results from 
these alternate mock catalogs, in addition to our fiducial one. 

As discussed in Sec. 2.1, the samples we use for our analysis 
are not volume-limited. We carried out our NFW fitting procedure 
to signals calculated using only 2 < 0.08 LBGs to quantitatively 
check what happens to the best-fitting masses. We found that with 
the blue 2 < 0.08 LBG sample, there was a slight tendency for the 
NFW masses to decrease by 15 per cent compared to our main re¬ 
sults, while for the red 2 < 0.08 LBG sample, the masses increase 
by typically 15 per cent. Thus, use of a volume-limited sample is 
not driving our conclusions about mass differences; indeed, with 
a 2 < 0.08 sample we would have found 30 per cent larger mass 
differences, albeit with substantially larger statistical errors as well. 

Another source of systematic uncertainty is the assumed cos¬ 
mology, which largely affects the values of concentration used in 
the fitting. We will consider this problem under the more general 
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Figure 11. Top: The mean halo masses inferred from NFW fits to AE 
as a function of the average stellar mass for each bin, shown separately 
for red (red circles) and blue (blue triangles) LBGs, along with shaded re¬ 
gions showing the 16-84 percentile range. The dashed lines show results 
without the mock catalog-based corrections. The single red and blue points 
shown separately with errorbars are the results for the single bin with M* > 
10^ ^Mq, which has a very broad stellar mass distribution that is different 
for the two color samples. Bottom: The difference in mass between the red 
and blue samples, expressed as logio{M 2 oom,red> - logio(A^ 200 m,blije>- 
The dark and light shaded regions show the 16-84 percentile and 2.5-97.5 
percentile ranges, respectively. 

umbrella of concentration-related uncertainties, which also include 
any errors in the adopted c{M) (even for the correct cosmology) 
and the impact of assembly bias (which could cause opposing er¬ 
rors in the concentrations for blue and red galaxies: red ones should 
reside in older halos which are more concentrated). In principle, 
some errors in the c{M) relation will he removed by the mock 
catalog-based correction factors. For example, if the c{M) relation 
used for the fits is wrong compared to both reality and the mocks, 
then our final results will be unaffected after applying the mock 
catalog-hased corrections. Since our primary conclusion is that red 
and blue central galaxies reside in halos with different masses on 
average, we focus particularly on the following question: could the 
impact of assembly bias on the concentrations of halos hosting red 
and blue central galaxies be causing us to incorrectly draw our pri¬ 
mary conclusion, rather than it being a real effect? 

If red galaxies actually live in halos with higher concentra¬ 


tions than we assumed and than in the simulations, and vice versa 
for blue galaxies, then the inferred halo masses for red (blue) galax¬ 
ies are too high (low). The sign of this effect is therefore correct to 
explain our result, but it is unclear whether the magnitude could be 
as large as the observed halo mass differences (factors of 2-2.5). 
To address this question, we note that Wechsler et al. (2002) found 
a scatter in logj^Q c of ~ 0.12. Taking the extreme assumption that 
red and blue central galaxies at fixed halo mass represent halo sam¬ 
ples split into the upper and lower part of the c(M) distribution 
depending on the red galaxy fraction, (cred) and (cbiue) can differ 
by up to ~ 35 per cent (depending on the mass range and red frac¬ 
tion, considering reasonable values of both quantities). Redoing the 
fits with modifications in concentration that are of that order for the 
stellar mass bin from 10^°'^ ^ M*/Mq < 10^^ causes the mass 
ratio to change by about 12 per cent. Even with extreme assump¬ 
tions, therefore, assembly bias cannot be responsible for any size¬ 
able part of the observed factor of 2-2.5 difference in mass between 
red and blue galaxies®. Under less extreme assumptions, there may 
be mass uncertainties at the few per cent level due to assembly bias 
affecting the halo concentrations for red and blue central galaxies. 
This is subdominant to the overall modeling uncertainties discussed 
above. 

As a final test of concentration-related uncertainties, we reran 
the fits with the concentration free but with a lognormal prior about 
the adopted c(M) relation used for the rest of this work. In gen¬ 
eral, the concentrations went down, and the masses went up. For 
red LBGs, we find that the masses typically increase by ~ 20 per 
cent in the data, and ~ 11 per cent in the mocks. Thus, when we 
consider the change in the mass correction factor that is based on 
the fits to the mock catalogs, our estimated red central galaxy halo 
masses would go up by typically ~ 9 per cent (but are also more 
noisy due to the extra degree of freedom). For blue LBGs, we find 
that the masses typically increase by ~ 30 per cent in the data, 
and ~ 20 per cent in the mocks. Thus, the net increase in the es¬ 
timated blue central galaxy halo masses when taking into account 
the changed correction factor would be ~ 10 per cent. Thus, our 
estimates for both blue and red central halo masses would increase 
by ~ 10 per cent, leaving the mass ratio unchanged on average, but 
noisier. The increase in noise and the lack of change in our con¬ 
clusion about halo mass bimodality justifies our choice to keep the 
fixed-concentration results as our primary ones. 

5.4 Comparison with previous results 

In this section, we compare our results with previous estimates of 
the galaxy type-dependent relationship between stellar and halo 
mass of central galaxies. In all cases when comparing with pa¬ 
pers that have multiple redshift ranges, we use the closest available 
range. Stellar masses have been converted to our adopted units and 
value of h and IMF, while halo masses have been converted to our 
adopted convention, M 2 oom with our value of Qm- We have not 
accounted for more subtle differences in halo vs. stellar mass rela¬ 
tions that could arise from, e.g., different assumptions about halo 
concentration vs. mass relations. Differences in stellar mass mod¬ 
eling can lead to additional (non-IMF-related) systematic offsets 
between stellar mass estimates that can be significant (Leauthaud 

® Also, for red galaxies, increasing the concentration by 35 per cent results 
in notably worse values; the data in fact prefer reduced concentrations 
rather than increased ones. For blue galaxies, there is less discriminating 
power between concentrations due to the lower signal to noise. 
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et al. 2012b), but are difficult to estimate and correct. When com¬ 
paring with previous SDSS results, we use the MPA/JHU stellar 
mass estimates that were adopted by many of those works. The ex¬ 
ception to this is More et al. (2011), who used stellar masses from 
the Bell & de Jong (2001) approach; however, the fitting formulae 
in Li & White (2009) show that for the range of stellar masses used 
in this work, the MPA/JHU and Bell & de Jong (2001) agree on 
average, so we do not apply a correction. In all cases, the plotted 
quantity is {M 2 oom) at fixed stellar mass, either in bins in stellar 
mass (shown using points) or from a parametric model or contin¬ 
uous function based on finely binned halos in simulations (shown 
using curves). The errorbars indicate error on the mean, rather than 
intrinsic scatter about the mean. 

A number of previous works measured the color-dependent 
average halo mass for central galaxies in bins in stellar mass. Man¬ 
delbaum et al. (2006), Velander et al. (2014), and Hudson et al. 
(2015) measured galaxy-galaxy lensing by galaxies split into pas¬ 
sive and active categories as a function of stellar mass, and used 
a simplified halo model to interpret the results in terms of typi¬ 
cal central halo masses (along with a satellite term). The top left 
panel of Fig. 12 shows a comparison between our results (shown 
as shaded regions that are the same as those in Fig. 11) and those 
three works. As shown, the results of Mandelbaum et al. (2006) are 
largely consistent with ours, with a hint of slightly lower masses 
for red galaxies below = 10^^ Mq. However, this could 

be a consequence of modeling differences (since a simple halo 
model was used to separate out satellite and central contributions, 
which are quite important for low-mass red galaxies). Alterna¬ 
tively, this is the only work against which we compare that used 
a morphological- rather than color-based separator (indicated on 
the plot using “early/late” instead of “red/blue”); while morphol¬ 
ogy and color are correlated, they do not split up the data in exactly 
the same way and could drive some of the observed differences. 
Also note that our current results have a slightly higher S/N than 
the earlier Mandelbaum et al. (2006) results, due to the larger area 
used and the better photometric redshifts for the source galaxies 
(which make our weighting closer to optimal). This S/N differ¬ 
ence, combined with the slight shift for the early-type galaxies, is 
the reason that we now see a more statistically significant difference 
between the average halo masses for red and blue central galaxies 
even below M* = lO^’^ Mq. The Velander et al. (2014) results 
seem consistent with ours for red galaxies, but somewhat high for 
blue galaxies. The Hudson et al. (2015) results are slightly lower 
than ours for red galaxies, but agree for blue galaxies. Differences 
in modeling and/or in stellar mass estimates may be at least partly 
responsible for these differences, which in any case are not very 
statistically significant. 

The lower left panel of Fig. 12 shows a comparison with two 
rather different results. More et al. (2011) used satellite kinemat¬ 
ics of central galaxies to infer the relation between halo and stel¬ 
lar mass. Their results agree with ours at high mass, but are sig¬ 
nificantly higher for M*'^*’^* < 10^^ Mq (however, at the low¬ 
est masses the discrepancy is not significant, primarily due to our 
much larger errors). One possible explanation for a discrepancy of 
this sign is satellite contamination of their central galaxy sample 
(Skibba et al. 2011; Lin et al. 2015), which could be more impor¬ 
tant below M* ~ lO^^M©; however. More et al. (2011) applied 
additional isolation criteria to avoid this problem. There are slight 
differences in their sample selection due to their different isolation 
criterion and color cuts, though it is unclear why those would cause 
a difference that always goes in the same direction. The most likely 
remaining explanations are modeling issues that may be more rele¬ 


vant at lower mass, where the number of satellites per halo becomes 
quite small on average; we defer further exploration to future work. 
For a comparison of several kinematics-based stellar vs. halo mass 
relations, see Wojtak & Mamon (2013). 

Also in the lower panel of Fig. 12 are the results of Rodrlguez- 
Puebla et al. (2015), who used a halo occupation distribution 
(HOD) to interpret the stellar mass function and galaxy cluster¬ 
ing of low-redshift galaxies in the SDSS. We show their mean halo 
mass as a function of stellar mass without errorbars, since the er¬ 
rors are very small compared to ours. Their results agree fairly 
well with ours for blue galaxies, but are below our results for red 
galaxies above 10^^ Mq. It is possible that in the absence of lens¬ 
ing, their halo masses are more reliant on assumptions within their 
halo model to interpret the clustering. Wang et al. (2015) show that 
even with the same cosmological model and stellar mass function, 
the clustering signal can vary significantly depending on model as¬ 
sumptions. We propose that this may be the primary cause for the 
discrepancy. 

The top right panel of Fig. 12 includes two papers that inter¬ 
preted galaxy-galaxy lensing and clustering measurements jointly. 
The results of Zu & Mandelbaum (2015a) were obtained by ap¬ 
plying a halo model to the galaxy-galaxy lensing and clustering 
measurements of Main sample galaxies divided by color and stel¬ 
lar mass (i.e., data similar to the “Main sample” curves in Figs. 8 
and 9). As shown, these agree extremely well with our LBG re¬ 
sults, despite the significantly more complex modeling needed to 
interpret those results. The halo model used there must account for 
the contributions by those Main sample galaxies that are satellites, 
and includes the effects of scatter between stellar and halo mass. 
The fact that our results are consistent with theirs after using our 
mock-catalog-based correction for that scatter in the LBG signals 
is a non-trivial result. 

The curves from Tinker et al. (2013) in the top right panel of 
Fig. 12 were obtained by jointly modeling the stellar mass function, 
lensing, and clustering in the COSMOS survey using a formalism 
that essentially doubles the number of halo model parameters in 
order to separately describe the results for red and blue galaxies. 
As shown, for blue galaxies their results are consistently slightly 
above ours, but with significant overlap of the error regions. For red 
galaxies, their results are consistent with ours for M, < 10^^ Mq, 
but deviate to significantly lower halo masses at high stellar mass. 
As a consequence of these two differences, the red vs. blue halo 
mass difference inferred from their results is smaller than ours. At 
the higher mass end, it is likely that their results are dominated by 
the higher signal-to-noise stellar mass function and clustering con¬ 
straints, which may be more dependent on modeling assumptions 
as mentioned above in the discussion of the Rodrlguez-Puebla et al. 
(2015) results. It is also worth noting that the mean redshift for this 
result is 0.36, so there may be some small evolutionary effects com¬ 
plicating the comparison between their results and ours. 

Finally, the bottom right panel of Fig. 12 shows the results 
of two mock catalogs. The first are the age-matching (Hearin & 
Watson 2013) mock catalogs publicly released with Hearin et al. 
(2014), who used stellar mass rather than luminosity as in the orig¬ 
inal age-matching paper. We show the mean halo mass at fixed stel¬ 
lar mass without errorbars, since the error on the mean is quite 
small. As shown, the {M 2 oom} curve for red central galaxies is 
slightly steeper than in our LBG results. However, as a basic con¬ 
sequence of their age-matching procedure, their {M 2 oom) for blue 
central galaxies are above the red ones for much of the stellar mass 
range. This results in a highly significant discrepancy with our ob¬ 
servations of blue central galaxy halo masses. The inversion of the 
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red and blue curves over a wide range of stellar mass, which results 
from the assumptions of the method and some particular implemen¬ 
tation details in Hearin et al. (2014), is discussed in more detail in 
Zu & Mandelbaum (2015a). It is worth noting, however, that Hearin 
et al. (2014) explored their ability to predict observable quantities 
associated with mixed galaxy samples containing central and satel¬ 
lite galaxies. In principle, a cancellation between effects seen in the 
two samples could be removed when exploring a sample of purely 
central galaxies, giving a more severe discrepancy with observa¬ 
tions in our work compared to the moderate tension shown in the 
plots in Hearin et al. (2014). As mentioned in Sec. 2.1, the fact 
that our sample is not volume-limited (and the one used by Hearin 
et al. 2014 is volume-limited) cannot be responsible for these differ¬ 
ences. A direct comparison of the AE profiles for all Main sample 
galaxies, and for LBGs, against predictions from the Hearin et al. 
(2014) is shown in Appendix C. 

The points shown in the bottom right panel of Fig. 12 are 
from the central galaxies in the semi-analytic mock catalogs used 
throughout this work (Sec. 3), using the same set of stellar mass 
bins as for the data. The statistical errors on the mean masses are 
not shown because they are very small; the systematic uncertainties 
are larger (Sec. 5.3). As shown, the curves for red central galaxies 
in the SAM-based catalogs are slightly above our observations, but 
not in a significant way except for at M* > Mq, also re¬ 

flected in Fig. 10. The deviation between our observations and the 
blue galaxy halo masses in the mocks is also more significant at 
the highest stellar masses (> 10^^ ATq) where we have very few 
galaxies. Overall, the mock catalogs exhibit a difference between 
the average halo masses for red and blue central galaxies in these 
bins that is consistently of the right sign and slightly smaller mag¬ 
nitude compared to the real data. This is remarkable since no clus¬ 
tering or lensing information was used when tuning the parameters 
controlling formation and quenching processes in this simulation. 

From a theoretical standpoint, it is interesting to consider how 
our results can be used. For example, the upper right panel of 
Fig. 12 shows a comparison against a halo model that includes a 
prescription for quenching (the shutdown of star formation activity 
that transforms a galaxy from active to passive) from Zu & Man¬ 
delbaum (2015a) and matches our results very well. This quench¬ 
ing model, necessary to make a color-dependent HOD, adds four 
parameters to the halo model, which is much fewer than the most 
obvious way to make a color-dependent HOD (which would dou¬ 
ble the number of parameters). Different models for the quenching 
process - i.e., different dependencies on stellar and/or halo mass - 
modify the mean halo mass at fixed stellar mass, suggesting that 
one can use our central galaxy constraints from this paper to con¬ 
strain quenching models. In contrast, the poor match to the average 
halo masses in the age-matching mocks from Hearin et al. (2014) 
calls into question the assumptions underlying the galaxy-dark mat¬ 
ter connection in those mock catalogs. 


6 DISCUSSION AND CONCLUSIONS 

In this work, we have used the galaxy-galaxy lensing signals of 
a sample of LBGs in the SDSS Main galaxy sample to measure 
the mean dark matter halo masses for central galaxies as a func¬ 
tion of their stellar mass and color. To validate the relatively sim¬ 
ple interpretation of these lensing signals in terms of average halo 
masses, we have used both simulation-based mock catalogs that 
closely reproduce the observed properties of Main sample galaxies 
and LBGs, and cross-checks within the observational dataset itself. 


We observe a clear, > 3a difference in the average halo masses 
of blue and red central galaxies at fixed stellar mass, for all stel¬ 
lar masses above 6 x 10^° ATq. The halo masses of the red (pas¬ 
sive) central galaxies are higher than those of blue (active) central 
galaxies. We have demonstrated that selection effects cannot be re¬ 
sponsible for this result, using a semi-analytic mock catalog that 
reproduces many properties of the LBG sample. 

Our result that there is a bimodality in halo masses of central 
galaxies at fixed stellar mass can be related to the underlying physi¬ 
cal processes such as feedback. For example, Wang & White (2012) 
used indirect measures of halo mass (the satellite counts around 
central galaxies) to infer a difference in average halo mass at fixed 
stellar mass for isolated red and blue primary galaxies. Their find¬ 
ing in a mock catalog that exhibited the similar effect (Guo et al. 
2011) was that this could be attributed to AGN feedback suppress¬ 
ing star formation in central galaxies that have massive black holes, 
thus lowering the average stellar mass for red galaxies at fixed halo 
mass (or, at fixed stellar mass, giving a larger average halo mass for 
red primaries than for blue ones). 

As discussed in Sec. 5.4, our results are largely consistent with 
previous lensing results (Mandelbaum et al. 2006; Velander et al. 
2014; Hudson et al. 2015), with only minor tensions that may re¬ 
sult from different modeling processes, stellar mass estimates, or 
ways of dividing the sample into blue and red subsets. We observe 
a clear tension with the results from satellite kinematics (More et al. 
2011) below M* = Mq, likely due to modeling differences. 
Our tension goes in the same direction for both the red and blue 
samples, with the kinematics results giving a higher mass; thus, 
our work and theirs come to the same basic conclusion that red 
galaxy samples have higher average halo masses than blue galaxy 
samples at fixed stellar masses. Our results for blue galaxies are 
consistent with those of Rodriguez-Puebla et al. (2015), who used 
an HOD model to interpret galaxy clustering and the stellar mass 
function, but our results for red galaxies are in disagreement with 
theirs, likely due to the increased sensitivity to modeling assump¬ 
tions when using clustering to infer halo masses. However, our re¬ 
sults are consistent with those based on HOD modeling of galaxy- 
galaxy lensing and clustering by Zu & Mandelbaum (2015a) using 
a formalism that incorporates the results of halo quenching. We 
also observe a highly significant disagreement with the blue cen¬ 
tral galaxy halo masses in the age-matching mocks from Hearin 
et al. (2014), due to an inversion in those mocks of the average halo 
masses for red and blue central galaxies. 

Our results clearly imply that high-fidelity mock catalogs that 
reproduce the galaxy-matter connection cannot be produced using 
methods that neglect type-dependence, such as simple halo model¬ 
ing. In addition, the age-matching approach appears to give average 
halo masses for blue central galaxies that are highly discrepant with 
our observations, as discussed in more depth in Sec. 5.4. More so¬ 
phisticated methods that include a strong galaxy type-dependence 
in {M 2 oom} (enforcing the correct sign of the difference) at fixed 
stellar mass seem to be necessary to reproduce observations; as we 
have demonstrated, this can be achieved with semi-analytic models 
or a type-dependent halo model. 
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Figure 12. All panels show the mean central halo mass as a function of stellar mass for red and blue galaxies from this paper (Fig. 11) as shaded regions. 
The top left panel shows results from measurements of galaxy-galaxy lensing by Mandelbaum et al. (2006), Velander et al. (2014), and Hudson et al. (2015); 
the bottom left panel shows results from More et al. (2011) (satellite kinematics) and Rodn'guez-Puebla et al. (2015) (HOD interpretation of abundance and 
clustering); the top right panel shows results from HOD modeling of lensing and clustering measurements by Zu & Mandelbaum (2015a) and Tinker et al. 
(2013); and the bottom right panel shows age-matching mocks from Hearin et al. (2014) as well as the semi-analytic catalogs used throughout this work.. All 
errors are la. 
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APPENDIX A: TESTS OF THE NON-VOLUME-LIMITED 
NATURE OF THE SAMPLE 

As described in Sec. 2.1, we compared the lensing signals for the 
full LEG samples (as shown in Fig. 7) with those for the LBGs at 
2 < 0.08. While technically the redshift below which the samples 
are volume limited depends on the stellar mass range, this serves 
as an approximate test for how much it matters that the samples 
are not volume-limited. Fig. Al compares the lensing signals for 
the red (top) and blue (bottom) LBGs with those for z < 0.08 
LBGs. The results are fully consistent with each other, supporting 
our choice to proceed with the full LBG samples. 


APPENDIX B: MASS RESULTS 

In Table B1 , we summarize our average halo mass results and other 
relevant statistics for red (top section) and blue (bottom section) 
central galaxies. The errorbars on the mean halo masses are actually 
the 16th and 84th percentile values of the non-Gaussian error dis¬ 
tribution. We have used those same percentile values to give non- 
Gaussian 16th and 84th percentile values for the ratio of stellar-to- 
halo mass (normalized by flb/Qm). Finally, the total corrections 
applied to the NFW masses, including the effects of completeness 
and impurity of the LBG sample as well as systematics in the fitting 
(e.g. due to the finite width of the halo mass distribution) are shown 
as the last column. These are the product of the two correction fac¬ 
tors derived in Sec. 3. 


APPENDIX C: AGE-MATCHING COMPARISON IN 
DETAIL 

In this appendix, we present the results of comparing AE profiles 
in the 2 = 0 catalog provided by Hearin et al. (2014) with the ob- 
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Figure Al. Comparison between lensing signals for all LBGs and for 2 < 0.08 LBGs for the red (top) and blue (bottom) subsamples, split into the same 
stellai* mass bins as in the rest of this paper. The data for all LBGs is the same as in Fig. 7. In both panels, the circles (triangles) show the results for LBGs 
(2 < 0.08 LBGs). Neither panel has 2 < 0.08 LBG results for the highest stellar mass sample, just because there are almost no lenses in that sample below 
that redshift. 
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Central galaxy halo masses 19 


Table Bl. Summary of our results for the average halo mass for central galaxies in fixed stellar mass bins for red (top section) and blue (bottom section) central 
galaxies. The quantities that are tabulated are the lensing-weighted stellar mass of the galaxies taking into account their weight in the lensing measurement 
for our canonical (VAGC) stellar masses as well as for the MPA/JHU stellar masses, the mean halo mass, the ratio of central stellar mass to host halo 
mass normalized by Q.i,/Q,rn assuming a Planck Collaboration (2014) cosmology for the VAGC and MPA/JHU stellar masses respectively, and finally the 
multiplicative correction factor used to estimate the mean halo mass for central galaxies from the best-fitting NFW mass. 


Intr In- / rff ^ , f (M 200 m)'\ M.,cen/M200. 
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10.86 

10.97 
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U.iUU_Q 

0 129+0 014 
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1.09 
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^•^^/_0.076 
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1.29 
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served AE profiles for red and blue LBGs and Main sample galax¬ 
ies in this work. This comparison avoids the mass inference proce¬ 
dure, and depends simply on replicating the LBG selection in those 
mock catalogs. While we use the z — 0 simulation snapshot, we 
weight the LBGs to reproduce the lensing-weighted stellar mass 
distribution in the observed (flux-limited) LBG sample. 

The results of the AE comparison are shown in Fig. CL Since 
our samples are not at 2 = 0, we must make a choice as to whether 
we take the AE profiles in their 2 = 0 simulation as being in 
physical separations, so no conversion is done to account for the 
differences in redshift; or whether their predictions are in comoving 
coordinates, requiring a conversion to compare with our results. We 
show both options in this figure, but the difference between them is 
not sufficient to substantially change our conclusions. 

As shown, for red galaxies, the age-matching predictions are 
only mildly discrepant with the data. At high mass, the mocks pre¬ 
dict significantly higher concentrations for red galaxies than are 
suggested by the observed AE profiles in the real data. The dis¬ 
crepancies between the model and data are particularly evident for 
the higher mass blue galaxies (middle and right panels), similar to 
the results shown in Fig. 12 when comparing masses. It is interest¬ 
ing to note that for M* < IC^Mq, the blue masses are higher 
than the red ones in Fig. 12, but this is not strongly reflected in the 
AE profiles for LBGs. This is primarily due to the fact that the 
model also includes a segregation of the red and blue into halos of 
different concentration, resulting in profile differences that partly 
counteract the effect of the mass differences and result in the red 
profiles being above the blue for < 0.2/i“^Mpc (and the oppo¬ 
site is true for larger Vp). For M* > Fig. 12 shows that 

the red galaxies have higher halo masses than the blue ones, and 
this is reflected in the AE profiles; however, the AE profiles for 
the blue LBGs (or all galaxies) are well above the data in this bin. 
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Figure Cl. Points with eiTorbars show lensing signals for red and blue LBGs (top) and Main sample galaxies without LBG selection (bottom) in stellar mass 
bins. The lines show the predictions from the Hearin et al. (2014) age-matching mock catalogs, with two lines for each sample (depending on whether we treat 
their predictions as being fixed in physical or comoving coordinates, with the latter always being higher than the former). 
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